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Hamilton's  principle  and  the  equations  of  motion  of  an  elastic  shell  with  and 
without  fluid  loading 
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Abstract 

It  has  proven  quite  difficult  to  employ  exact  elastodynamic  theory  to  describe  the  behavior  of 
elastic  vibrations  on  arbitrary  bounded  shells.  In  addition,  exact  theories  preclude  direct 
interpretation  of  particular  features  observed  due  to  the  excitation  of  elastic  shell  surfaces.  A 
rather  interesting  approach  to  describe  surface  vibrations  may  be  obtained  by  constructing  a 
Hamiltonian  in  some  approximate  form  that  assumes  some  correlation  of  motion  of  the  outer 
and  inner  shell  surface.  The  class  of  theories  that  allow  for  this  approach  are  referred  to  in 
applied  mechanics  as  shell  theories.  The  interesting  feature  of  this  Hamiltonian  approach  is  that 
one  can  add  various  physical  mechanisms  to  the  Hamiltonian  such  as  extensional  motion,  rotary 
inertia,  fluid  loading,  etc.,  and  thereby  study  the  individual  contributions  to  resonance  patterns 
while  adding  physical  insight  to  the  fundamental  processes  that  occur  on  shell  surfaces.  We 
develop  shell  theories  in  this  manner  and  examine  various  contributions  via  Hamilton's 
principle.  We  believe  that  fluid  loading  has  by  and  large  not  been  treated  adequately  in  the  past, 
and  we  place  particular  emphasis  on  the  treatment  of  that  contribution  to  this  work. 


1.  INTRODUCTION 

Sound  scattering  from  submerged  elastic  shells  is  of  interest  to  a  broad  community  of 
scientists  and  engineers.  For  example  in  the  area  of  nondestructive  testing  ultrasonic  scattering 
is  of  considerable  interest  to  the  structural  engineering  community,  while  in  the  area  of  remote 
sensing  or  object  identification  scattering  from  objects  proves  invaluable.  The  acoustical 
scattering  from  elastic  objects  or  more  simply  the  generation  of  resonances  from  bounded  elastic 
shells  involves  the  use  of  acoustic  and  elastodynamic  theory.  Exact  solutions  for  scattering  from 
elastic  targets  exist  for  target  shapes  for  which  the  elastodynamic  equation  is  separable.  For 
three  dimensional  targets  separability  is  only  possible  for  spherical  and  rectangular  targets.  For 
more  complicated  targets  numerical  techniques  must  be  used.  For  axi symmetric  smooth  targets 
the  extended  boundary  method  of  Waterman  has  proven  useful  though  limited  and  an 
approximate  theory  would  be  desirable  for  more  complicated  targets.  The  purpose  of  this  work 
is  to  develop  an  approximate  theory  based  on  shell  theory  that  would  be  of  use  for  general 
shapes  that  could  also  include  structural  loading.  Shell  theories  afford  a  very  powerful 
methodology  to  "build  in"  structural  features  usually  via  a  variational  principle  in  which  some 
Largrangian  or  Hamiltonian  is  constructed  by  introducing  physical  features.  The  simplest 
theories  which  always  assume  correlated  motion  of  the  inner  and  outer  surfaces  of  the  target 
include  kinematic  features  and  potential  energy  terms  based  on  the  generalized  Hooke's  law 
which  leads  to  the  lowest  order  symmetric  mode.  One  can  add  rotational  inertia  to  the 
Hamiltonian  which  allows  an  antisymmetric  mode.  However  the  antisymmatric  mode  with  just 
this  added  feature  does  not  have  die  correct  asymptotic  behavior.  Timoshenko  noted  that  by 
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including  shear  distention  the  antisymmetric  mode  obeyed  the  correct  asymptotic  features.  One 
goes  on  from  there  to  include  higher  order  corrections.  Further,  proper  fluid  loading  —  which 
from  our  experience  is  no  easy  task  --  must  be  introduced  correctly.  The  proper  inclusion  of 
fluid  loading  introduces  a  fluid  borne  wave:  the  pseudo-Stoneley  wave.  Clearly  thin  shell 
theories  are  rather  geometrical  and  must  be  constructed  for  each  shape  due  to  their  dependence 
on  shape  dependent  dynamic  factors.  Non-spherical  three  dimensional  objects  present  a 
problem  which  we  seek  to  address  in  the  future.  In  this  work  we  wish  to  develop  and  test  a 
suitable  shell  theory  for  spherical  shells  that  can  be  readily  generalized  to  spheroids  and  finite 
cylinders.  In  the  following  section  we  outline  the  derivation  of  a  shell  theory  with  many  of  the 
featured  described.  The  theories,  with  various  levels  of  sophistication,  are  then  used  to  compare 
with  the  exact  method.  The  results  are  then  discussed  and  future  work  is  described. 

2.  DERIVATION  OF  EQUATIONS  OF  MOTION 


In  spherical  shells  membrane  stresses  (proportional  to  (})  predominate  over  flexural  stresses 
(proportional  to  (3 2)  where 


/*  = 


1  h 

Vl2  a 


(1) 


We  differ  from  the  standard  derivation  for  the  sphere  [2]  by  retaining  all  terms  of  order  /J2in 
both  the  kinetic  and  potential  energy  parts  of  the  Lagrangian  and  by  considering  the  resonance 
frequencies  for  the  fluid  loaded  case  to  be  complex.  We  note  that  this  level  of  approximation 
will  allow  us  to  include  the  effects  of  rotary  inertia  in  our  shell  theory,  as  well  as  damping  by 

fluid  loading.  The  parameter  /?  itself  is  proportional  to  the  radius  of  gyration  of  a  differential 
element  of  the  shell  and  arises  from  integration  through  the  thickness  of  the  shell  in  a  radial 

direction.  We  will  use  an  implicit  harmonic  time  variation  of  the  form  exp(-ia)t).  We  begin 
our  derivation  by  considering  a  u,v,w  axis  system  on  the  middle  surface  of  a  spherical  shell  of 

radius  a  (measured  to  mid-shell)  with  thickness  h,  where  u  increases  meridionally  southward, 
v  increases  latitudinally  eastward,  and  w  increases  radially  outward. 

2.1.  Lagrangian  Variational  Analysis 
Our  Lagrangian,  L,  is 

L  =  T-V  +  W,  (2) 

where  T  is  the  kinetic  energy,  V  is  the  potential  energy,  and  W  is  the  work  due  to  the  pressure 
at  the  surface.  The  kinetic  energy  is  given  by 

T  =  KIC  +  *,2  Xfl  +  *)2  sin  Odxdddtp,  ( 3 ) 


where  the  surface  displacements  are  taken  to  be  linear: 


,,  xs.  x  dw 
u,  =  (l  +  — )u — — , 
a  a  dd 


(4) 
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and 


wt  =  w. 


(5) 


The  motion  of  the  spherical  shell  is  axisymmetric  since  the  sound  field  is  torsionless.  Thus 
there  is  no  motion  in  the  v-dircction.  Substitution  of  Eqs.  (4)  and  (5)  into  Eq.  (3)  yields,  after 

integration  over  x  and  <p, 

_  f*  ■  nr,  k*  h 1  2  2  /j5  h*  dw 

T  =  npA  sin#[( - -  + —  +  ha  )u  -  2( - r  + — )u— - 

80a2  2  80a2  ^  2f> 
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+( 


h5  . h 3 


■  +  ~)(^-y+(:r-  +  ha2)w2]dd. 


80 a2  12  '  <90"  12 


(6) 


or,  in  terms  of  P, 


T  =  np,ha2 JjO.804  +  6p2  +  l)zi2  -  (3.6>34  +  6 p2)u 


+(1.8  p*  +  p2){^-)2  +(J32  +  l)vv2]sin  Odd, 
d6 


.  dw 

le 


(7) 


where  the  first  and  last  terms  in  square  brackets  in  Eq.  (7)  are  associated  with  linear 
translational  kinetic  energies  and  the  middle  two  terms  are  associated  with  rotational  kinetic 
energies  of  an  element  of  the  shell. 

The  potential  energy  of  the  shell  is 

1  f*  f2*  fA/2  - 

V  =  2  Jo  Jo  J-A/2(0w£w  +  +  0)  ^  ddXCidd<t>'  <  8  ) 


where  the  nonvanishing  components  of  the  strain  are 


r*+j+±i 

rdu 

<?V 

\  89  )  a 2 

JO 

~  do2, 

(9) 


and 


1 


eM  =  —(cot  8u  +  w)  + 
a 


and  where  the  nonzero  stress  components  are 
E 


(e„  +  veH), 


(10) 


(11) 
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~  J  ^  ^  £ H  + 


where  E  is  Young's  modulus.  By  substitution  the  potential  energy  becomes 

,,  1  f* f2* f*/2  T  E  1  ( r/i  x\^u  xd2w  2  x.  xdwA  l2 

V  =  —  - r - r[(l  +  -)- - — T-  +  w2  ]  +{cOt  0[(1  +  — )U - — +  w) 

2  Jo  Jo  j-k/2  l-v2(x  +  a)2v  a  dd  add 2  a  add 


.  .  x.  xdw.  ....  x.du  x  d2w  }  ,2  . 

+2 v{cot 0[(1  +  —  )u - ]  +  w}[(l  +  — )— - —r +  w]  (x  +  a)  smddxdd,  (13) 

a  a  dd  a  dd  a  dd  J J 

which  after  integration  is 

_  J}Eh  r  +  +  (w  +  u  t  +  2  y(w  +  — )(w  +  ucot  d) 

l-v2**  £0  00 

o2 r,du  d2w. 2  2  dw~  _  dw..du  d2w„.  .  ..  .. 

+p  km~W}  cot  +^c«9(U--)(--^r)]|smede.  (14) 

Terms  in  the  potential  energy  proportional  to  /J2  are  due  to  bending  stresses. 

And  finally,  the  work  done  by  the  pressure  of  the  surrounding  fluid  on  the  spherical  shell  is 
given  by 


IV  =  2na2  f  p.wsinddd, 

JO 


where  pm  is  the  pressure  at  the  surface. 

2.2.  The  Lagrangian  density  and  its  equations  of  motion 

Integration  along  the  polar  angle  d  is  intrinsic  to  the  problem,  therefore  we  must  turn  to  a 
Lagrangian  density  formulation  to  solve  for  the  equations  of  motion.  Our  Lagrangian  density  is 
just 

L  =  *p>2[(l  +  602  +  l&F)u2-(6p2  +  3.6p4)u~+(P2+l.ty*)A2 

dtt  du 

+(1  +/J2)w2]sin  d  -  {(w  +  -^i)2  +  (w  +  «cot  d)2  +  2  v(w  +  —  )(w  +  ucot  d) 

1  —  v  dd  dd 

02r ,du  d2w.2  2  n,  2  ^  dw..du  d2w...  . 

+l>  +CM  +2wco,e(u--)(--^)i)sm# 

+2«a2p.wsin  0, 


(16) 
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with  corresponding  differential  equations 


0  = 


dL  d  dL  d  dL 
du  dd  du9  dt  du,  ' 


(17) 


and 


.  dL  d  dL  d  dL  d2  dL  d2  dL 
0  =  — - — - —  + - - - +  ■ 


dw  d6  dwg  dt  dw,  dddt  dwg,  dd2  dwM ' 

where  subscripts  denote  differentiation  of  the  variable  with  respect  to  the  subscript 
By  substitution  of  Eqs.  (17)  and  (18)  into  (16)  we  obtain 


(18) 


0  =  (1  +  fi2)\  ~~ + cot  0 -  ( v + cot2  0)u 
A  du 


+t(l  +  V)  +  02(  v+cot2  0 )]0  -  0(1. 80*  +  6p2  + 1)0  -  (1.8 P4  + 


(19) 


and 


d*w  „  n,  „  d3w 


+cot  0[(2  -  v + cot2  0)f}2  -  (1  +  v)]u  -  pi2  -  2f}2  cot  0 

d2w 


d04  d0* 


+0 2(1+  v+cot  0)~^j~0' cot  0(2-  v  +  cot  0)—  -2(1+  v)w 


t1  o)~ 

dd1^  r  ”  d0 

+0-(i.8/34  +  ~  d-8)34  +  3/32)cot 00 


(20) 


These  differential  equations  of  motion  (19)  and  (20)  have  solutions  of  the  form 


»0 


(21) 


and 


>v(r?)  =  X^(»7). 


•■0 


(22) 
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where  rj  =  cos 9  and  Pa(t})  are  the  Legendre  polynomials  of  the  first  kind  of  order  n.  When 
the  differential  equations  of  motion  (19)  and  (20)  are  expanded  in  terms  of  Eqs.  (21)  and  (22), 

we  obtain  a  set  of  linear  equations  in  terms  of  Ua  and  W„,  whose  determinant  must  vanish.  We 
shall  consider  two  cases:  with  and  without  fluid  loading. 

2.3.  Vacuum  case 

The  simpler  case  is  that  when  the  spherical  shell  is  surrounded  by  a  vacuum  such  that  there  is 

no  damping.  In  this  case,  the  pressure  at  the  surface  vanishes:  pt  =  0.  The  set  of  linear 
equations  the  expansion  coefficients  must  satisfy  are 

0  =  [Q2(l  +  6/J2  + 1.8/)4)  -  (1  +  P2)k]U„  +  [Q2(3/?2  + 1.804)  -  p2K- (1  +  v)]W„,  (23) 

and 

0  =  -A.[( k  -  3 )p2  - 1.8/34  + 1  +  v]l/.  +  [Q2(l  +  2 p2  + 1.804 ) - 2(1  +  v) -p2KX JVV.,  (24) 

where  Q  =  (oa/  cp,  k  =  v+ X„  - 1,  and  A„  =  n(n  + 1).  In  order  for  Eqs.  (23)  and  (24)  to  be 
satisfied  simultaneously  with  a  non-trivial  solution  the  determinant  of  the  system  must  vanish: 

0  =  a4(l  +  6 p2  + 1.804)(1  +  2  p2  + 1.8  pA) 

+a2{(3 p2  +  l.SpA)X,[(K-3)p2  -l.&p*  +  l+  v] 

42(1  +  v) +p2  vA.](l  +  6P2  + 1.804)  -  (1  +  p2)K{\  +  ip2  + 1.8/34 )} 

+(1  +  P2 ) k[ 2(1  +  v)  +  p2 kX„ ]-X„[(k-  3 )p2  -  l.Sp4  + 1  +  v](Jp2K  + 1  +  v).  (25) 

Since  there  are  no  damping  terms,  the  shell  vibrates  theoretically  forever.  Thus,  the  normalized 
frequency  Q  can  be  taken  to  be  real.  Equation  (25)  is  quadratic  in  Q2,  thus  we  expect  two  real 
roots  to  (25)  and  thus  two  modes  for  the  motion  of  the  shell.  They  are  the  symmetric  and 
antisymmetric  modes. 

2.4.  Fluid  loaded  case 

For  the  fluid  loaded  case,  we  must  consider  a  modal  expansion  of  the  surface  pressure  in 
terms  of  the  specific  acoustic  impedance  z,.  In  its  most  general  form  this  is 


••  M 

p(a,9,<p)  =  IX*  X./r(cos  6)cosm<f>, 

M>0  «a0 


(26) 


where 


.  hSka) 
z.  =  ipc~ — 

K(ka) 


(27) 


The  specific  acoustic  impedance  za  can  be  split  into  real  and  imaginary  parts: 
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zm=rm-i(om„  (28) 

where 


r„  =  pc  Re 


ih,(ka) } 

Km  j  ’ 


and 


(29) 


m. 


ih'(ka)} 

Km’' 


(30) 


For  the  case  we  are  considering  of  axisymmetric  motion,  the  surface  pressure  is  given  by 


p.<0>= cose), 

<i*0 


(31) 


or 


P.(8)  =  -  (o2WmrnK)Pmm6).  (32) 

»*0 

Use  of  Eq.  (32)  in  our  set  of  differential  equations  of  motion  (19)  and  (20)  yields  the  following 
set  of  linear  equations  for  the  expansion  coefficients  in  the  case  of  a  fluid  loaded  spherical  shell: 

0  =  [Q2(l  +  6/J2  +  1.8/34) -(1  +  /J2  )*•]£/„  +[Q2(302  +  1.804)-J32kt-(1+  v)]Wa,  (33) 

and 

0  =  -Aj(<r-  3)f}2  -1.804  +  1+  v]Um 

+{Q2(1  +  a  +  202  +  1.8/J4)-  2(1  +  v)  +  Qiy -&2kX  JWm,  (34) 

where 


a  = 


(35) 


and 
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£Jk_ 

*  P.c, ' 


(36) 


Again  the  determinant  of  Eqs.  (33)  and  (34)  must  vanish.  However,  in  this  instance  the  value 
of  Q  must  be  taken  to  be  complex;  the  resonances  have  a  width  that  depends  on  the  damping. 
The  result  of  setting  this  determinant  to  zero  is 


0  =  Q4(l  +  6/J2  + 1.804)(1  +  a  +  2p2  + 1.804) 

+Q3/y(l  +  6j82  +  1.804) 

+Q2  {Op2  + 1. 804 ) A„[(  k  -  3)p2  - 1. 8/J4  + 1  +  vl 

-{2(1  +  v)  +  /?2#trA.](l  +  6p2  + 1. 8/34 )  -  (1  +  p2 )  v(l  +  a  +  202  + 1. 8  p* )} 

+n[-iy(l  +  p2)K] 

H\+p2) x{2(l  +v)  +  P2kX„\-X'[(k- 3 )p2  - 1. 8/34  + 1  +  v)(p2 k  + 1  +  v).  ( 37 ) 


Equation  (37)  has  four  complex  roots.  From  work  with  an  exact  modal  solution  to  the 
problem,  we  expect  two  roots  to  be  associated  with  the  symmetric  and  antisymmetric  modes  of 
the  shell.  We  expect  the  other  two  roots  to  be  associated  with  a  water-borne  pseudo-Stoneley 
wave. 

3.  CONCLUSIONS 

The  next  step  is  to  plot  the  roots  of  Eqs.  (25)  and  (37)  to  compare  the  resonances  predicted 
by  these  models  with  those  given  by  exact  modal  expansion  solutions.  By  suppressing  a  and 
y,  the  model  associated  with  Eq.  (37)  reverts  to  the  vacuum  case  model  associated  with  Eq. 

(25).  Similarly  suppression  of  factors  of  P  in  Eq.  (25)  will  result  in  a  reversion  to  a  previously 
derived  solution  [1].  We  may  then  rank  the  three  different  models  according  to  their  degree  of 
physicality  and  compare  their  results  for  various  relative  shell  thicknesses  against  each  other  and 
against  the  exact  results  of  the  modal  expansion  model.  We  may  also  consider  the  limitations  of 
each  of  the  models  including  the  exact  solution,  as  well  as  those  of  shell  models  in  general. 

These  models,  fluid  loaded,  vacuo  case,  and  membrane,  are  successively  less  physically 
sophisticated  and  give  successively  less  good  comparison  with  exact  (modal  expansion)  results. 
Starting  with  the  least  sophisticated  model,  we  see  in  Fig.  1  thick  spherical  steel  shell 
dilatational  (symmetric)  and  flexural  (antisymmetric)  mode  resonances  calculated  by  the 
membrane  model.  Here  and  in  the  succeed:ng  figures  thick  means  hla  =  0.1;  thin  means  hta  = 

0.01.  The  shell  material  is  a  generic  steel  with  density  p,  =  7.7  times  that  of  water,  shear 

velocity v,  =  3.24  km/s,  and  longitudinal  velocity  v,  =  5.95  km/s.  The  surrounding  fluid  is 

taken  to  be  water  with  density  p  =  1000  kg/m^  and  sound  velocity  cw  =  1.4825  km/s.  The 
symmetric  mode  shows  a  good  comparison  between  exact  and  shell  theory  predictions,  but  the 
antisymmetric  shell  theory  results  for  this  approximation  compare  poorly  with  the  exact  flexural 
results.  Note  that  some  symmetric  mode  resonances  were  not  found  by  our  exact  theory 
algorithm.  In  Fig.  2  we  see  thin  spherical  steel  shell  dilatational  (symmetric)  and  flexural 
(antisymmetric)  mode  resonances  calculated  by  the  membrane  model.  Again  there  is  good 
comparison  between  dilatational  (symmetric)  mode  resonances  calculated  by  the  two  methods, 
except  for  the  first  couple  of  resonances.  The  exact  flexural  resonances  are  highly  damped;  the 
few  resonances  seen  are  water-borne  pseudo-Stonely  resonances.  And  again  the  shell  theory 
flexural  (antisymmetric)  mode  resonances  do  not  asymptote  properly  with  increasing  order.  In 
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Fig.  3.  we  have  thick  spherical  steel  shell  dilatadonal  (symmetric)  and  flexural  (antisymmetric) 
mode  resonances  calculated  by  shell  theory  without  fluid  loading  (vacuum).  As  in  the 
membrane  model  the  shell  theory  and  exact  calculations  compare  well  for  the  dilatadonal 
(symmetric)  mode  resonances.  In  contrast  with  the  membrane  model,  however,  the  exact  aid 
shell  theory  calculations  for  this  model  show  much  better  agreement  for  the  flexural 
(antisymmetric)  mode  resonances.  This  model  does  not  include  fluid  loading,  but  does  include 
the  effects  of  rotary  inertia.  The  vacuum  shell  theory  flexural  mode  resonances  do  not 
asympto  for  large  size  parameter  ka  to  the  exact  results,  however.  In  Fig.  4.  we  see  thin 
spherical  steel  shell  dilatational  (symmetric)  and  flexural  (antisymmetric)  mode  resonances 
calculated  by  shell  theory  without  fluid  loading  (vacuum).  As  in  the  membrane  model  the  shell 
theory  and  exact  calculations  compare  well  for  die  dilatational  (symmetric)  mode  resonances 
except  for  the  first  couple  -  >f  resonances.  This  vacuum  model  does  not  have  fluid  loading,  and 
has  insufficient  damping  for  the  first  two  dilatational  (symmetric)  mode  resonances.  Again,  the 
flexural  (symmetric)  mode  resonances  show  roughly  the  correct  behavior,  but  it  is  not  possible 
to  tell  what  the  asymptotic  value  of  the  phase  velocity  would  be  for  large  size  parameter  on  this 
scale.  Next  in  Fig.  5.  we  have  a  plot  of  thick  spherical  steel  shell  dilatational  (symmetric)  and 
flexural  (antisymmetric)  mode  resonances  calculated  by  shell  theory  with  fluid  loading.  As  in 
the  vacuum  case  as  well  as  for  the  membrane  model,  the  dilatational  (symmetric)  mode 
resonances  compare  well  for  exact  and  shell  theory  methods.  The  flexural  (antisymmetric) 
mode  resonances,  as  calculated  by  shell  theory  with  fluid  loading,  do  not  have  the  correct 
asymptotic  limit  for  large  size  parameter,  although  they  do  exhibit  the  correct  behavior  for  lower 
values  of  ka.  Finally,  in  Fig.  6.  we  see  thin  spherical  steel  shell  dilatational  (symmetric)  and 
flexural  (antisymmetric)  mode  resonances  calculated  by  shell  theory  with  fluid  loading.  The 
exact  and  shell  theory  calculations  agree  well  for  the  dilatational  (symmetric)  resonances  and 
exhibit  a  marked  improvement  for  the  first  several  shell  theory  symmetric  mode  resonances. 
This  is  due  to  the  inclusion  of  fluid  loading  in  the  model.  The  flexural  (antisymmetric)  mode 
resonances  show  the  appropriate  behavior  on  this  rather  limited  size  parameter  scale. 

The  proper  correction  for  the  large  ka  asymptotic  resonance  behavior  is  most  easily  found  by 
including  the  shear  stress  distortion  along  with  a  Timoshenko-Mindlin[2,3]  shape  factor. 
Inclusion  of  the  shear  distortion  in  the  potential  energy  would  make  the  flexural  modes 
asymptote  to  the  coincidence  velocity  and  the  shape  factor  can  be  adjusted  so  that  they 
asymptote  to  the  Rayleigh  velocity  as  expected. 
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